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Abstract 



A significant challenge in zero-field /uSR experiments arises from the uncertainty in the muon site. It is possible to 
calculate the dipole field (and hence precession frequency v) at any particular site given the magnetic moment pt and 
magnetic structure. One can also evaluate /(v), the probability distribution function of v assuming that the muon site 
can be anywhere within the unit cell with equal probability, excluding physically forbidden sites. Since v is obtained 
from experiment, what we would like to know is g(fi\v), the probability density function of [i given the observed v. 
This can be obtained from our calculated f(v/fj) using Bayes' theorem. We describe an approach to this problem 
which we have used to extract information about real systems including a low-moment osmate compound, a family of 
molecular magnets, and an iron-arsenide compound. 
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1. Introduction 

In a fiSR experiment the muon-spin precession frequency v can be used to deduce the local magnetic field B at 
the muon site. If there are a number of muon sites with different local fields {Bj}, then the juSR signal can contain a 
number of components with frequecies {vy}. In complex systems it can be highly non-trivial to determine the precise 
location of the muon site or sites. However, it is nevertheless useful to attempt to extract an estimate of the magnitude 
of the moment of the magnetic species producing the local field, even in the face of ignorance of the location of the 
muon site or sites. In this paper we describe a method to attempt this using Bayesian inference. The paper is structured 
as follows: in Sections 2 and 3 we review the theory of dipolar fields and Bayes' theorem respectively to provide the 
necessary background to the calculation which is described in Section 4. 

2. Dipolar fields 

An implanted muon spin precesses around a local magnetic field, Bi oca i, with a frequency v = (y /J /2^)|Bi oca i|, 
where y fi /2n = 135.5 MHz T -1 . The magnetic field Bi oca i at the muon site is given by 



^local - #o + ^dipole + B L + B t 
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where Bo represents the applied field (zero in the experiments considered here), Bdip is the dipolar field from magnetic 
ions, Bl = f*oM/3 is the Lorentz field, Bdemag is the demagnetizing field from the sample surface and Bhyperfme is the 
contact hyperfine field caused by any spin density overlapping with the muon wavefunction. In antiferromagnets the 
Lorentz and demagnetizing fields vanish (in polycrystalline ferromagnets they cancel to some extent). The contact 
hyperfine field is hard to estimate but we will neglect it. The remaining term is the dipolar field Bd; p and is a function 



of the muon-site r^. It can be written as 



SdipM-Z^V,)™?, (2) 



a sum over the magnetic ions; the magnetic moment of the ith ion is m,. In Eq. D°f(r^) is the dipolar tensor given 
by 

r 3RfR^ 



R z 
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(3) 



where B, = (Rf, R y r Bf) = r jl - r, and &# is the Kronecker delta (6 a/s = 1 if a = p, else = 0). The behaviour of this 
tensor is dominated by the arrangement of the nearest-neighbour magnetic ions and leads to a non-zero local magnetic 
field for almost all possible muon sites, even in an antiferromagnetically ordered system fl]|2]. The sum in Eq. Q is 
taken over the infinite lattice, but it is well known [3| that this sum converges in such a way that it is necessary only 
to sum over points inside a sphere centred on r u with sufficiently large radius. An alternative method of calculation is 
provided by the method of Ewald summation (for details see [4]). 



3. Bayes' theorem 

We recap some elementary probability theory J5]|6l. The conditional probability P(A|B) is the probability that 
event A occurs given that event B has happened. The joint probability P{A n B) is the probability that event A and 
event B both occur. The joint probability P(A n B) is equal to the probability that event B occurred multiplied by the 
probability that A occurred, given that B did, i.e., 

P(A n B) = P(A|B)P(B), (4) 

and, equally well, 

P(A n B) = P(B|A)P(A). (5) 
Now consider the case where there are a number of mutually exclusive events A, such that 

Y J P(A i ) = l. (6) 

Then we can write the probability of some other event X as 

P(X) = Y j P(X\A i )P(A i ). (7) 

i 

In very general terms, one can say that given some hypothesis H there usually exists some computational strategy to 
evaluate the probability of a particular outcome O assuming that hypothesis to be correct (i.e., there is some method 
to compute the quantity P(0|H)). However, what you often want to do is the reverse of this: you know the outcome 
because it has actually occurred and you want to choose an explanation out of the possible hypotheses. In other words, 
given the outcome you want to know the probability that the hypothesis is true, and the problem is that P(H|0) is 
typically much more challenging to evaluate. The needed transformation of /"(OIH) into P(H|0) can be accomplished 
using Bayes' theorem (named after Thomas Bayes (1702-1761), although the modern form is due to Laplace). This 
theorem can be stated as follows: 

W(B|A) 

Here P(A) is called the prior probability, since it is the probability of A occurring without any knowledge as to the 
outcome of B. The quantity which you derive is P(A|B), the posterior probability. The proof of Bayes' theorem is 
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very simple: one simply equates Eqs. Q and |5]) and rearranges. For the purpose of this paper, we will write Bayes' 
theorem using Eq. fT) as 

mpm 

P(P\V) = -r. ■ (9) 

jP{v\y.')P(n')dn' 



4. Calculation 

Since the positive muon seeks out areas of negative charge density, constraints can be placed on the likely location 
of stopped muons. For example, the muon is unlikely to stop close to the positively charged ions in a system. In 
many oxides, muons have been shown to stop around 0.1 nm from an 2 ~ ion |7|. In our calculations, we assume a 
magnetic moment fi on the magnetic species in our material and consider a particular magnetic structure. Positions 
in the unit cell are then generated at random and, provided the relevant constraints are satisfied, the dipole field 
is calculated at each of them. After many such randomly generated positions, one obtains a distribution of dipole 
fields. (This distribution sometimes have sharp features associated with them which are van Hove singularities 018].) 
The magnitudes of the resulting dipole fields are then converted into muon precession frequencies, and the resulting 
histogram yields the probability density function (pdf) f(v/fi), evaluated as a function of precession frequency v 
divided by magnetic moment /u (since the precession frequency scales with the magnetic moment). This function 
/(v/yu) allows us to evaluate P(v\fi) of Eq. |9]). We can write 

P(v\fi) = -f(v/»). (10) 

The function f(v//j) is normalized so that J f(v/fi) d(v//j) — 1, and hence the factor of ^ is needed in Eq. |Io| ) so 

that f P(v\fi)dv = 1. An example of this approach is shown in Fig. 1(a) for a cubic lattice of antiferromagnetically 
coupled spins. The solid line shows the case when a simple constraint is applied so that the muon is not permitted to 
stop at a site closer than a particular critical distance to the magnetic moments, resulting in a cut-off of the tail at high 
frequency. 




Figure 1: (a) The dipolar field distribution for a simple cubic lattice of antiferromagnetically aligned magnetic moments parallel to [110] (with 
antiferromagnetic wave vector q = (n, n, 7r)). The lattice parameter is a and the size of the moment is /i |2 |. The solid (dashed) line shows the 
distribution with a cut-off (not) applied, (b) The extracted pdf for the moment given a particular observed frequency (v = 2nfift/a 3 ). 
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Since v is obtained from a real experiment, what we would like to know is g(/-i\v), the pdf of yt given the observed 
v. This can be obtained from our calculated f(v/fj.) using Bayes' theorem in the form of Eq. |9]i, which yields 

jT £/(v//i')<fc' 

where we have assumed a prior probability [P(/u)] for the magnetic moment that is uniform between zero and // max , 
and so P(/S) is replaced by the uniform probability density 1 /yu max [which cancels on the top and bottom of Eq. (jTTJ] . 
We choose /z max to take a large value, although we have found that our results are insensitive to the precise value of 
A<max- A very simple example of this approach is shown in Fig. 1(b). When multiple frequencies v, are present in the 
spectra, it is necessary to multiply their probabilities of observation in order to obtain the chance of their simultaneous 
observation, so we evaluate g(/i|{v,}) oc ]~[, f(Vi/(J) dv,, where Av,- is the error on the fitted frequency. 

We have now applied this technique to /iSR data a variety of real systems in which the muon site is not known. 
These include Ba2NaOs06 in which we can show from the observed precession frequencies that the magnetic ground 
state is most likely to be low-moment (« 0.2 fj.%) ferromagnetism and not canted antiferromagnetism |9|. We have 
also used it to show a reduced moment in the two-dimensional molecular magnet [Cu(HF2)(pyz)2]BF4 [ 10 1 and in the 
pnictide superconductor NaFeAs Hill . In all these cases we do not have a priori information concerning the muon 
site but can nevertheless place bounds upon the magnetic moment from the observed precession signal using this 
technique. A possible drawback that should be borne in mind is that the hyperfine contribution to the local field is 
neglected and if this is significant it could affect the conclusions drawn. As many of the systems examined so far 
using this technique have localized, reduced moments and lower-frequency precession signals, it is probable that the 
hyperfine contribution is not significant in these cases. 
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